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Abstract 

The perturbation expansion for a general class of many-fermion systems with a non-nested, 
non-spherical Fermi surface is renormalized to all orders. In the limit as the infrared cutoff 
is removed, the counterterms converge to a finite limit which is differentiable in the band 
structure. The map from the renormalized to the bare band structure is shown to be locally 
injective. A new classification of graphs as overlapping or non-overlapping is given, and 
improved power counting bounds are derived from it. They imply that the only subgraphs 
that can generate r factorials in the r^^ order of the renormalized perturbation series are 
indeed the ladder graphs and thus give a precise sense to the statement that 'ladders are 
the most divergent diagrams'. Our results apply directly to the Hubbard model at any 
filling except for half-filling. The half-filled Hubbard model is treated in another place. 
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1. Introduction and Overview 



1.1 The Problem 

Consider the following problem in many-body physics. Let A be a finite box in d- 
dimensional space, i.e. A C IR^ or A C F, where F is a lattice in IR^, and let Co-(x) 
and c^(x) be fermionic annihilation and creation operators obeying the canonical anti- 
commutation relations {co-(x), cjr,(x')} = 5o-ct'<^(x — x') and let be the fermionic Fock 
space generated by this algebra [BR] . Let H\ = Ho + XV be the operator on given by 

Ho = J2 f c?s(x)c+(x)(T + C/)c,(x) (LI) 

where T is an operator describing the one-particle kinetic energy, U is multiplication by 
a periodic potential, and J ds{x) denotes Jj^^dx for a continuous system and X^^eA ^ 
system on a lattice. Let no-(x) = c+(x)ccr(x) be the number operator at x for spin a. The 
interaction 

V = J2j J ds{^')n^{x.)v^^,{^ - x')n<,.(x') (L2) 

a,a' 

is assumed to be short-ranged (see Assumption Al below). The Hamiltonian H\ describes 
many electrons in a crystal or on a lattice, that interact with a stationary ionic background 
through U and with each other through the pair potential V. If the coupling strength of 
the electron-electron interaction A = 0, the electrons move independently according to 
the one-particle Schrodinger operator T + U{x). In the continuum system T = —A/2m 
is the Laplacean and U{x + 7) = U{x) for all 7 G F, where the lattice F is generated 
by d linearly independent vectors in IR'^ (e.g. F = 2^^); in the case of a lattice system, 
U = and the kinetic energy T is defined by the hopping matrix between the sites of 
the lattice. For A 7^ 0, the potential V takes into account interactions such as screened 
electromagnetic interactions. A slight generalization of (1.2) allows for inclusion of phonon- 
mediated interactions. 
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Let j3 = 1/kT be the inverse temperature and define the grand canonical partition 
function Z\ as 

Za= tr e-^^'^^-'^^^) (1.3) 

where 

Na = J2 f ds{^)na{^) (1-4) 

- A 

is the number operator, /j, is the chemical potential and the trace is over Fock space. For 
an observable O, i.e. a polynomial in the fermion operators, the thermal expectation value 
is defined as 

{0)a = 7^ tr f e-'^^"'^-^'^^^0) (1.5) 
The question wc are interested in is whether the thermodynamic limit Q = lim Q\ of the 

A^oo 

connected Green functions Qa = (c^^ (x^) . . . Ca'^ (x^)) a, conn , which are special cases of O 
above, exists and whether in infinite volume a weak-coupling expansion 

oo 

g = J2yGr (1.6) 

r=0 

can be used to determine the dependence of Q on A. 

For this question the most interesting, because most singular, case is that of zero tem- 
perature, T = 0. For positive temperature or the finite volume lattice case the expansion 
obtained by expanding the factor e^^ in A is convergent, but its radius of convergence 
shrinks to zero in the thermodynamic and zero-temperature limit: at T = and in infinite 
volume, one can not even pose the question of convergence of the expansion in A because 
the coefficients Gr already diverge for r > 3. In the limit T ^ 0, (1.5) reduces to expecta- 
tion values in the ground state of the system, so physically the question is about the nature 
of the many-particle ground state of the system and the validity of perturbation theory to 
calculate n-point- functions. The radius of convergence of the unrenormalized expansion in 
finite volume shrinks to zero as the volume goes to infinity. Thus, although the expansion 
converges for the large but finite systems which these models are to describe, this is true 
only if A is of order 1/volume, which is obviously unrealistic for any macroscopic system. 
Consequently, the unrenormalized expansion will not give insight into the properties of the 
ground state. 

In this paper we consider formal perturbation theory. That is, we study the ther- 
modynamic limit of the coefficient functions Gr- By an analysis similar to [FTl], the 
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expansion is renormalized so that these functions converge as the volume goes to infinity. 
More precisely, we introduce a well-defined infinite volume model obtained by cutting off 
the singularity at the Fermi surface (i.e. introducing an infrared cutoff) and renormalized 
by including counterterms K in the action, and then show that all coefficients Gr have 
limits as the infrared cutoff is removed. These counterterms are bilinear in the fermions 
and can therefore be viewed as a modification of Hq (although they are treated as extra 
interaction vertices in the formal expansion). They also have finite limits as the infrared 
cutoff is removed. The limiting counterterms reflect the modiflcation of the band structure 
due to the interaction. The precise meaning of this will be discussed in much more detail 
below. Although we do not go through the finite-volume bounds here, it will be clear 
from the way our bounds are derived that the same procedure can be applied to obtain an 
expansion in finite volume with coefficients that converge in the thermodynamic limit. 

Except for special cases, the renormalized expansion is, as an expansion in A, not 
convergent but only locally Borel summable because the coefficients behave as Gr ~ r\. 
The occurence of these factorials indicates that the nonperturbative ground state may 
exhibit symmetry breaking. For example, if the interaction is attractive in the zero angular 
momentum sector, this is the case [FT2]. One of the main results we shall prove here is 
that for a very wide class of models (and regardless of the sign of the interaction), the r 
factorials in individual graphs come only from ladder diagrams. 

Renormalization has been done in [FTl] for the continuum case where T = — A/2m 
and t/ = 0. We shall refer to this case as the spherical case since the band structure (defined 
below) has an 0{d) rotational symmetry. The procedure for removing the divergences in 
the present case is similar to the spherical case in that we have to renormalize two-legged 
insertions. However, the present work is a nontrivial extension of [FTl] because in contrast 
to the spherical case the counterterms are not constants. In brief, subtracting functions is 
much more complicated than subtracting constants. In particular, the regularity properties 
of the counterterms are quite subtle. 

In the remainder of this introductory section, we give a non-rigorous, physical discus- 
sion of why divergences occur and how they may be removed by renormalization. We hope 
that this will convince the reader, before going through all the details, that the renormal- 
ization subtractions are natural and the divergences of the naive expansion are artificial 
in these models. We state our main results in Section 1.5 and then discuss their physical 
interpretation. Finally, we give an overview of the sections containing the proofs. Every 
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section begins with a brief explanation of what is done and how it fits into the general 
strategy. 



1.2 The Formal Perturbation Expansion 

The models have the formal functional integral representation 

P(r;, ^) = y" L)V^V;e^+(^'^)+('''^) (1.7) 

where A= — (V', C~^'ip) — XV, DipDip is the formal measure Ylx a dipa{x)d'ipa{x), 

(V5,C-V) = / ds{x)ds{y)J2M^) {C-')^^{x,y)My). (1-8) 

a, 13 

and 

V= ds{x)ds{x') ^ ■4^a{x)'^l3{x)Val3,a',(i'{x,x')-4)a'{x')'^(i'{x'), (1.9) 

where now / ds{x)F{x) stands for the integral over the spatial variable x and imaginary- 
time r, x= (t, x), with an appropriate measure, e.g. 

ds{x) = drd'^x (1.10) 

for a continuous system on [0, (3] x R*^ and 

I ds{x)F{x)= f dT^F{T,:i^) (1.11) 

J Jo 



for a lattice system on [0, 0] x F, e.g. F = Z*^. Here /? = l/ksT is the inverse temperature. 
The imaginary time is introduced to get a functional integral representation for the trace 
over Fock space in the standard way. The connected Green functions can formally be 
calculated as derivatives of log P with respect to the sources r] and fj. 

In this paper, we consider the limiting case T = 0, so /? = oo and the configuration 
spaces are JR'^'^^ and H x F (e.g. H x Z*^), respectively. The spin index a e {T)i}) the 
interaction is assumed to be translation invariant, so that 

Val3,a'p' {x, x') = Va/3,a'/3' (t - t', X - x'), (1.12) 
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and short-range, i.e. v decreasing so fast that its Fourier transform v is at least twice 
differentiable (see Assumption Al below). Note that we do not assume that it is instan- 
taneous. For simplicity, we also assume that it is spin-diagonal, i.e. Va(3,a'(3' = ^aa'Sp/S'V- 
In contrast to the assumption about the decay of v, the latter assumption is merely for 
notational convenience and can easily be dropped. 

One may imagine v to arise from exchange of (quasi) particles like photons or phonons 
and formalize this by a Hubbard-Stratonovich transformation, introducing one or more 
scalar fields with covariance v so that the interaction vertex is resolved as an exchange of 
fields and the interaction becomes bilinear in the fermion fields. For the purposes of the 
perturbation expansion we shall not need this. In particular since we assume smoothness 
of V, we shall not need a cutoff on the interaction lines, and we shall often draw graphs 
with four-legged vertices instead of ones with interaction lines. 

For the lattice models, we take 

(C-l) {X,X') = 6af3 ((^xx' {dr' - //) - Tx_x') S{t - t'), (1.13) 

where fi is the chemical potential and Tx_x' is the amplitude for hopping from site x to 
site x', which we assume to be symmetric and short- ranged (see Assumptions A2 and A3 
on e below). 

A model of particular interest that is easy to formulate but difficult to analyze is the 
Hubbard model, for which 

Tx = J] ty5x,y (1.14) 

|y|=i 

with ty the so-called hopping parameters. In the simplest version of the model, ty = t is 
the same for all y of length one, so the operator T is just the discrete Laplacean on 2^, 
with the diagonal term omitted since it can be absorbed in the chemical potential //, and 
the interaction term is on-site and spin-diagonal, 

Vaf3,a'f3'i^ - x') = 5xx'^a/3(^a'/3' • (1-15) 

Various extensions of this model, e.g. with more complicated finite range hopping have 
been studied in connection to high-temperature superconductivity. For suitable values of 
the filling factor, they all fall into the class of band structures discussed here. For a review 
of mathematically rigorous results about the Hubbard model, see [L] . 
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Formally equivalent to P, but in fact much more convenient is the generating func- 
tional for connected amputated Green functions 



where the constant Z takes out the field-independent term so that ^(0,0) = 0. as 
written above, is not a well-defined object in infinite volume; it can be made well-defined 
by restricting to a finite volume A, or by introducing a suitable cutoff. If the free covariance 
C is bounded and any power of it is integrable, j^^A exists and is analytic in A, as was first 
observed by Caianiello. However, for any realistic model, C will not have this properties, 
unless cutoffs are imposed. The radius of convergence obtained using naive bounds shrinks 
to zero when the cutoffs are removed, and establishing analyticity uniformly in the cutoffs 
requires techniques as in [FMRT]. 

Our analysis is done in momentum space, where from now on momentum is short for 
Bloch's quasi-momentum, which can be used to label one-particle states because of the 
periodicity of the one-particle potential U. In infinite volume, momentum space is the 
first Brillouin zone B, i.e. the torus 



where is the dual lattice to F, e.g. F* = 27r2 for F = 2 . In finite volume, the 
momenta are in a finite subset of i3, p = 27rn/L with n G 2*^ HB if the volume is a box of 
sidelength L. The eigenfunction expansions used to transform into momentum space are 
discussed briefiy in Appendix B for the general case; for the purposes of this introduction, 
we just give the formulas for the case of a lattice model on Z'^, where we can simply do a 
Fourier expansion. The only changes in the general case are (of course) that the Brillouin 
zone will differ with the lattice and that the formulas for switching between position and 
quasi-momentum space involve the eigenfunctions of the one-particle Hamiltonian Hq with 
the periodic potential. Under the Fourier transform 




(1.16) 



B = H'^/F* 



(1.17) 




(1.18) 



the quadratic part of the action becomes 




(1.19) 
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where we have dropped the hats and introduced the band structure 



e{p) = s{p)-fi (1.20) 

where 

£(p) = J ds(x)e-^P^T,, (1.21) 
and the interaction becomes, with pi = {ipi)o,Pi), and 

^-. = ...^p=^|^, (1.22) 



(1.23) 



V = J ^+Vl---^+V4(27r)^+l5((p2+P4-pi-p3)o)5#(p2-pi+P4-p3) 
^(P3 - Pl)V'(Pl)V'(P2)V'(P3)V'(P4). 

here 5^ is the delta function on B, more exphcitly 

'^''(P) = 72^ ^ E'^(P + 7) (1-24) 

where the 5 on the right side denotes that on H''. In general, the solution of the one-particle 
problem will produce crossing bands. We exclude this case here, and we also introduce 
an ultraviolet cutoff that removes the high energy bands. For the lattice systems, such 
a cutoff is already built in as the lattice spacing; for continuous systems it is not a real 
physical restriction since high energies do not occur in a crystal. If there arc finitely many 
bands that do not cross, the band index is just a bookkeeping device dragged along, so, 
without loss, we restrict to the one-band case here. 

For A = 0, the fermions do not influence each other and the model is completely 
characterized by the covariance C, 

/„-ipoT+ipx 
Wo - e(p) 

in the sense that all 2n-point functions are determinants of matrices with elements C {xi — 

Xj). 

The propagator in momentum space, C{p) = e'^°^^/(zpo — c(p)), has a singularity at 
Po = for all p E S, where /S = {p : e(p) = 0} is the Fermi surface of the independent 
electron approximation. Although the function l/{ipo — e(p)) is in lI^J{JR x B) for all 
(5 e [0, 1), graphs in the perturbation expansion diverge because of the singularity on S and 
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because in the expansion, arbitrary powers of C are integrated. The numerator e*^°°^ is 
included in the standard way since we want to consider the expansion around the situation 
where all states inside the Fermi surface, i.e. those with e(p) < 0, are already occupied. 
Expanding ^ in a formal power series in A, we can write 

6;(V'e,V5e) = 5]^'^r(V'e,V5e) (1-26) 
r>0 

with 

/2m / m 2m \ 

m^i r-'- i=l \i=l i=m+l / 

— a^,...,am 

m 

(G'2m,r)c,^,...,a^ (^'l' " " ■ ^ P2m-l)W'>Pai{Pi)i^ai{Pm+i) , 



(1.27) 



i=l 



where the coefficient function G2m,r is totally antisymmetric in the simultaneous exchange 
of momenta and spin indices (see Section 2.3). Again, the is periodic with respect to 

in the spatial part of the momentum. The coefficient G2m,r can be expressed in the 
usual way as a sum over values of connected Feynman diagrams. The sum over m runs 
over a finite index set for each fixed r because the number of vertices is r and the graphs 
are connected with 2m external legs. 

The Feynman graphs are similar to those in quantum electrodynamics: there are two 
types of lines, namely fermion lines (drawn solid), carrying a direction, and interaction lines 
(drawn dashed). The vertices have two legs to which fermion lines can be connected (one 
incoming, one outgoing), and one leg for an interaction line. The action determines the 
assignment of propagators C{p) to fermion lines, v{p) to interaction lines, and momentum 
conservation delta functions to vertices. Equivalently, one can replace two vertices that 
are joined by an interaction line by a single four-fermion vertex with exactly two incoming 
fermion legs and exactly two outgoing fermion legs. The graphs then have only four-lcggcd 
fermion vertices and only fermion lines. There is one notable difi^erence between the cases 
U = and C/ 7^ 0: In the spherical case (t/ = 0), where e{p) = pV2m, p G M'^. The 
corresponding ultraviolet problem (behaviour at large |p|) was solved in [FTl]. In presence 
of a crystal potential {U 7^ 0), the integrals over the spatial part of the momentum are 
over the first Brillouin zone B, which is a compact set. Thus there is no case of large p 
here. Momentum conservation at every vertex means conservation in B, as given by 
above. If one prefers to think of the momenta in H^, fixing momenta with 5^ means that 
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at every vertex, there remains a sum over 7 e F^. Although formally infinite, this sum 
always contains only one nonzero term since there is a unique 7 e that translates back 
a vector in WC^ into the fundamental domain of the translational group F'^. However, it is 
natural and simpler to consider momentum space as the torus B since e is F^-periodic. 
For example, in the Hubbard model, 

d 

e(p) = 2t ^ cos Pi - // (1.28) 

i=l 

is the tight-binding band relation and v{p) = 1. 

The much more general class of models and the range of chemical potential fx that we 
treat in this paper is given by the following assumptions. 



1.3 Assumptions 

We assume that the one-particle problem (discussed in Appendix B) is such that we have 
a Brillouin zone B which is a d-dimensional torus of type (1-17). We assume that e = e — ji 
(see (1.20)) is a continuous function on B and that for some value ji^ of the chemical 
potential, the Fermi surface 

S = {peB: e(p) = 0} (1.29) 

has only a finite number of connected components. Furthermore, there is A; > 2 and a 
neighbourhood N ol S such that: 

Al The interaction v G C^(1R x i3, (D). The sup norm over IR x i3 of the first k 
derivatives is finite. 

A2 The band structure e e C^{M, H), and Ve(p) for all TpeS. 

The third assumption is a geometrical condition on the Fermi surface. It is very simple to 
understand and is fulfilled for generic surfaces. Let n : S ^ H*^, u 1— > n(a;) = |^(<^), be 
the unit normal to the surface. By A2, 5' is a C'^ submanifold of B, and n is a C''~^ unit 
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vector field. If S consists of more than one connected component, choose a normal field 
for any component. For u!,uj' e S, define the angle between n(a;) and n{uj') by 

9{(jj,(jj') = arccos(n(a;) • n(a;')). (1.30) 

Let 

V(oj) = {u;' eS: \n{u;) ■ n(oj')\ = 1} = {u' E S : n{oj) = ±n(oj')}, (1.31) 

and denote the {d — l)-dimensional measure of A C 5" by vold-iA. Also, for any A C R'^ 
and /3 > denote by Ui3{A) = {p G : distance(p, A) < (3} the open /9-neighbourhood 
of A. For fixed e and Ho, we assume: 

A3 There is an open interval Ai around fio and there are strictly positive numbers 
Zq, Z^, p, l3o and k such that for all G the Fermi surface S = S{n) = {p G 
B : e(p) = 0} has the following properties: S{iJ,) C jV, and for all (3 < (3o and all 
a; G 5, 

{i) vold-i{Up{V{u;))nS)<Zof3- 

(ii) if a;' ^ Uf3{V{u;)) n <S, then |sin^(a;, a;')| = ^/l - {n{uj) ■ n{uj'))^ > Z^pp. 

Throughout this paper, A1-A3 will be assumed to hold, and will be assumed to lie in 
the interval M. specified in A3. We now explain what these assumptions mean. 

Assumption Al on {) is a decay assumption in position space, e.g. for an instantaneous 
interaction V on a lattice system on Z'^ and k — 2, Al holds if 

J2 |xnF(x)|<oo. (1.32) 

For continuous systems, Al is implied by a similar integral condition. 

Assumption A2 excludes singular points. For example, a point p on 5 where Ve(p) = 
is called a van Hove singularity. 

The condition that e is continuously differentiable is fulfilled for the case where e comes 
from a Schrodinger equation for the one-body problem with a regular periodic potential, 
if there is no level-crossing. Indeed, it is real analytic. In lattice models with finite-range 
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hopping, e is analytic. However, infinite-range hopping is also allowed: e e if the k*'^ 
moment of the hopping amplitude exists, i.e. J2 l^l^l'I'xl < oo. 

X 

Assumption A3 is, more informally, that for every uj & S 
(i) the set of points co' where the normal n{uj') is parallel or antiparallel to n(a;), has 

positive codimension k > in S and 
(m) if uj' is not in the set T>{lj), where the normal is (anti)parallel to n{u)), the angle 

between n{u>) and n{u>') increases with some power of the distance between u;' 

and V{u>). 

Thus in order to violate these assumptions, the surface S must have flat regions or 
subsets where 9{uj,u}') vanishes exponentially fast as \u} — uj'\ 0. To illustrate A3, we 
draw an example of a Fermi surface that satisfies A3 m. d = 2 (i.e. a Fermi curve) on 
B = ^llTx-n? (the square bounds the fundamental region [— tt, tt)^ for the torus and 
the shaded areas indicate e(p) < 0): 
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